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Abstract: We consider the effective theory for the low-energy modes of SU{2) 
'. gauge theory on the three-sphere. By explicitely integrating out the high- 

^ I energy modes, the one-loop correction to the hamiltonian for this problem 

Tj- I is obtained. We calculate the influence of this correction on the glueball 



spectrum. 



^ ■ 1 Introduction 

^1 Non-perturbative effects in gauge theories can be studied using finite volumes [jl], ||, |[. In 
■ a small volume asymptotic freedom implies that the coupling constant is small: this means 
that we can use standard perturbation theory. By increasing the volume, we can study the 
^ . onset of non-perturbative phenomena. The effects that we want to study are related to the 
^ ! multiple vacuum structure of the theory. For increasing volume, the wave functional starts 
to spread out over the configuration space in those directions where the potential energy 
is lowest, i.e. in the direction of the low-energy modes of the gauge field. In particular, it 
will flow over the instanton barrier that connects gauge copies of the vacuum. 

As long as the non-perturbative effects manifest themselves appreciably only in a small 
number of low-lying energy modes, this can be described adequately using a hamiltonian 
formulation. It is hence our strategy to split up the gauge field in orthogonal modes and to 
reduce the dynamics of this infinite number of degrees of freedom to a quantum mechanical 
problem with a finite number of modes. 

In our approach, we impose the Coulomb gauge by restricting the gauge fields to a 
so-called fundamental domain H. The spreading out over configuration space means 
that the wave functional will become sensitive to the boundary conditions that have to be 
imposed on the boundary of this fundamental domain. In these boundary conditions the 
dependence on the ^-angle will show up. 
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For more details on this method the reader is referred to where we used the lowest 
order effective hamiltonian in a variational calculation of the spectrum. In the present 
letter we perform the one-loop calculation, that is, we integrate out the high-energy modes 
in the path integral to obtain the correction to the lowest order hamiltonian. We will 
subsequently use this new hamiltonian to find the glueball spectrum. 

2 The effective theory 

We will briefly review the technical set up for the analysis on the three-sphere. For details 
and more motivation, we again refer to |0 and references therein. Let be the normal 
vector on the three-sphere. We define two orthonormal framings on by 

where we used the 't Hooft t] symbols J^. These 77 symbols occur in the multiplication 
rules for the unit quaternions: 

where and their conjugates = a^^ are defined by 

a^ = a,^r), a^ = (l,_zf). (3) 

We choose to write a gauge field on {n^A^ = 0) with respect to the framing e^: 

A, = Ae^ = A"e;^ (4) 

We introduce a number of su{2) angular momentum operators. We define = L\ = ^ ^^fi^^ 
and = I e^^d^. These operators generate the SU{2) x SU{2) symmetry of S"^ and satisfy 

L\ = L\. We introduce a spin operator 5* by {S°'A)i = —iSaijAj and an isospin operator T 
by T" = ad (Ta/2). We also define Ki = + Si and Ji = Ki + Ti. 
To isolate the lowest energy levels, we write 

^cM) = j,, \ = tr(A.-M,A,) + O {A) , (5) 

where the quadratic fluctuation operator M. can be rewritten as 

Ai« = (i?- - p);; (6) 

The zero-modes of M. correspond to pure-gauge modes of the gauge field. The 18 dimen- 
sional space B{c,d) given by 

B,{c, d) = (c«e; + d]e^) y = (c^ + ^1^/) (7) 

is the eigenspace of Ai corresponding to its lowest positive eigenvalue 4, whereas the next 
eigenvalue is 9. The tunnelling path is = —uS"', d^=0 with u running from to 2. For 
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ti = 1 it passes through the sphaleron, which is a saddle point of the energy functionaL 
The energy functional for these 18 modes is given by 

^cl(C' d) = I, \ = ^cl(c) + ^cl(^) + l(tK^) tr(y) - tr(Xy)), (8) 

Vcl(c) = 2tr(X) + 6detc+ i (tr^{X) - ir{X^)) , (9) 

with the symmetric matrices X and Y given by X = c and Y — d(f^. The lowest order 
hamiltonian for these modes is 

■2 

with / = and R the reinstated radius of the sphere. 



3 Gauge fixing 

We will impose the background gauge condition on the high-energy modes. Consider a 
general gauge field on 5"^: 

A^^{Ao,Ai), (11) 

where Aq is the time component of the gauge field and Ai are the space components with 
respect to the framing e^. Wc will now project out the background field B{c, d). Let Ps be 
the projector on the constant scalar modes, and let Py = Pc + Pd be the projector on the 
(c, (i)-space. We define the background field B and the quantum field Q by respectively 

= {PA)^^{PsAo,{PvA)i), (12) 
= A^-B^. (13) 

We define the gauge fixing function x by 

X = (1 - Ps) D^{PA)A^ + PsAo. (14) 

We use X fo impose the background gauge condition: x = is equivalent to Bq — and 
Dij^[B)Q^ = 0. We perform the standard manipulations with the partition function: after 
introducing Faddeev-Popov ghosts and expanding the classical action up to second order 
in Q we obtain 

Z = [ DBkD'Q^D'^D'i^exp \ ( U {-D^{B){1 - P)D^{B)}'il^ 
+ \fI^{B) - 2{D,F,,){B)Q, + Q,W,,iB)Q, 

with 

1^00 =-Dl{B) 



(15) 



Wo^ =-W^,o = -2ad(5,) (16) 
W,, = -2ad (F,,(B)) - {DliB)),, + 26,, 
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Remember that the covariant derivative Di{B) acting on vectors (or tensors) gives extra 
connection terms (due to 5*^ being a curved manifold), e.g. 

(A(5)C), = d,C, + [B,,C.j\-e,,kCu (17) 
= {-2iU + iBtT'' -iS,)^^Ck (18) 

The primed integration means that we have excluded the (c, (i)-modes from the integration 
over the vector field Qk-, and the constant modes from the integration over the scalar fields 
V', and Qo- 

The action contains a term J^Q^ with = {D^F^^){B). Since B need not satisfy 
the equations of motion, this term does not vanish. When expanding the path integral in 
Feynman diagrams, this term will give rise to extra diagrams, where J acts as a source. It 
can be shown that J will only contribute to terms in the effective lagrangian that we will 
consider to give only small corrections: they are at least of the order c^rf^ or c^d"^ and we 
will ignore them. 

Dropping the term linear in Q^, we obtain from Z the effective action 

= sijc^[B] - Indet' {-D^{B){1 - P)D,{B)) + \ Indet' {W^,) (19) 



4 The effective potential 

As the one-loop computation is the central result for this paper we present here a few of 
the details that are crucial to understand why we were able to complete this calculation. 
For computing the effective potential, B is considered to be independent of time. This 
effective potential is determined up to fourth order in the fields c and d and up to sixth 
order in the tunnelling parameter m, where we expect the physics to be most sensitive to 
the precise shape of the potential. We use A{B) to denote any of the operators occuring 
in ([T9|). They are of the form 

A{B) = -dl~dl + A{B). (20) 

In order to neatly perform the dimensional regularisation, we introduced a laplacian term 
for an ^-dimensional torus of size L that we attached to our space 0. The scale L should 
of course be chosen proportional to the radius R of the three-sphere. The assumption 
[(9o, A{B)\ = allows us to write 

where we took the time periodic with period T. The operator A can be expressed in terms 
of the angular momentum operators defined above, the functions V/ and the constants 
and d°;. To take the trace, we need a basis of functions. For the scalar operators (the ghost 
operator and Wqq), we can use \lmL m/j)(S)|l m^) or equivalently \lmR; j rrij), where m^, m/j 
and nij correspond to the z-components of Li, L2 and J respectively. Here / = 0, -1, 1, . . . 
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and j = \l — 1\, . . . ,1 + 1. For the vector operator Wij, we can use \lmR] knik) ® |1 rrit) or 
\lmji; k]j rrij), where the bounds on the various quantum numbers are obvious. Note that 
the c and d modes correspond to the vector modes with {I, k) = (0, 1) and (/, k) = (1, 0) 
respectively. For the scalar operators, the trace must not be taken over the / = functions, 
whereas for the vector operator the trace must not include the c and d modes. Note that for 
the case of the ghost operator, the operator (1 — P) (see eq. (p!9D) causes some intermediate 
vector modes to be projected to zero. 

Using the appropriate basis, we can calculate the determinants exactly for the vacuum 
u = (5 = 0), and for the sphaleron configuration u = 1 (c" = —5" and d = 0). The final 
summation over / is expressed using the function ({s, a), which is defined by 

Cis,a) = Y:^jj^^, Re(.)>i, (22) 

After analytic continuation we have 



oo 



as,a)=Y.^-^i-^riCRi2s + 2m)-l), .^1,-1,..., (23) 

m=0 



where denotes the Riemann ^-function and {s)m is Pochhammer's symbol. If s ap- 
proaches one of the poles, this expansion can be used to split off the divergent term: the 
remainder of the series is denoted with (f{s, a)- We find for the effective potentials 

V«(i? = 0) = -18 + 3Cfi(-3)-3Ci?,(-l), (24) 
V^^HSph) = -l-3v^-12v^+|v^-fv^+g-^7 + ^log(2) 

+T log (277) + -3) + CF(-f , -2) + Cf{-1 1) 

+3 CF(-i -3) + Cf(-|, -2) - 5 CF(-i 1) (25) 
The pole term is absorbed through the usual renormalisation of the coupling constant 



4e 

1 _ 1 11 
7r ~7o^ l^'e 



— = — + ^ ^_r,_ + finite renormalisation (26) 



For a general configuration along the tunneling path, we make an expansion of the eigen- 
values of the operators in u. Although it is still possible to calculate the spectra exactly, 
a polynomial form of the effective potential is much more useful in the variational method 
that we will use. We used Bloch perturbation theory to obtain the expansions of the 
eigenvalues of Wij. Results up to tenth order in u were obtained, where we used both 
MATHEMATICA and FORM in the calculations. 

As can be seen from fig. 1, the expansions do not converge to the exact result at m = 1. 
This should come as no big surprise, since we have no reason to expect the radius of 
convergence of the expansion to be as large as one. To find the effective potential for larger 
u, we write u = 1 + a and make a similar expansion of the effective potential around the 
sphaleron. Using the fourth order expansion in u and the first order expansion in a (i.e. 
the value and the slope of the potential at the sphaleron), we can construct a polynomial 
in u of degree six that is a good approximation to the effective potential. 
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We now turn to B = B{c) (we still keep d = 0). The perturbative evaluation of 
the individual eigenvalues is no longer possible, but we can use the following technique. 
Suppose we have A = F + A with F = A{0) such that [F, A] = 0. This allows us to 
substitute in (pi]) 

tr (e-^^) = tr (e-^^e"^^) = ^ e~^^^ J2 i^^) " ^^'^^ 



n=0 



The sum over z is a sum over the eigenspaces of F, Fi is the corresponding eigenvalue 
and trj () denotes a trace within the eigenspace. For both the scalar operators we have 
F = 4L^ and [L^, A{c)] = 0. The remaining problem of calculating tr; (^A^{c)^ then reduces 
to calculating traces of the form 

tTi{Li,...Li„T''' ...T""-) (28) 

which can be done relatively easy. For the vector operator we have F = 2L^ + 2K^ and 
[K'^, W{c)] 7^ 0. The trace in (^) can however still be written as a sum of traces in the 
different eigenspaces of F: 

tr (e-^^) 
= ^e-^^Hri(e-^(^-^'+^) 



Ee-^f:^tr.({(F-F.) + i}' 

E ^Pm (29) 



Here sp^ is the contribution of order m in A. Let Pi denote the projector on the eigenspace 
of Fi, and let T be given by 

Ti = tri(i) = tr (Pii), (30) 
T,j... = tiiPiAPjA---). (31) 



functions. With Aji = Fj — Fi we find 



For a given value of m we can perform the combinatorics to write sp^ in terms of the T 
= Fj — Fi we find 

sPi = Ee-^^'(-^)^^> (32) 

i 

E E tr, {A{F - F,r-'A) 



SP2 



„=2 



ij n=2 "'■ 



oo / \n n— 3 

SP3 = En-'^e-^'^^E^ E A,TAL-'-"^\ (34) 

ijk n=3 ' mi=0 
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etc. 

Starting from one eigenspace i^, the number of intermediate states that can be reached 
is finite. This allows us to extract one overal summation, and to perform the remaining 
finite sums. For the case at hand, wc extract the summation over I and are left with an 
/-dependent expression in which the T functions have the form 

tr,(P,^iyP,^iy---). (35) 

Since the intermediate modes in the T functions can only have A; = / — 1, ...,/ + !, we can 
write 

= ao{k) + ai{k)K^ + a2{k){K^)\ (36) 

for certain values of the coefficients. Note that a T function with an intermediate value 
A; = should be discarded, since this corresponds to a d mode. The combinatorics for sp^, 
the finite sums, the evaluation of the T functions and the final summation over I were all 
done in FORM. 

For general B{c,d), not even commutes with the various operators. The methods 
described above are however sufficient. The T traces that we have to calculate now also 

— * — 

contain the operators L2 and explicit V functions, as well as projectors on different L 
intermediate levels. Using the c d symmetry it is however only necessary to obtain the 
c^ii^ terms. Since the precise form of the coefficients is not very illuminating, we postpone 
writing down the effective potential until we have performed the renormalisation. 



5 The renormalisation 

To obtain the one-loop contribution to the operator B^, we perform the usual expansion of 
the path integral in Feynman diagrams. The subtleties related to the summation over the 
space-momenta were dealt with in the previous section. The diagrams needed are depicted 
in fig. 2, where the particle in the loop is respectively a ghost, a scalar Qq or a vector 
particle. The insertions in the diagrams correspond to the operator A defined above. The 
propagators in momentum space are given by 

(37) 



for scalar particles, and 



ki + k^, + + 1) 

kl + kl + 2/(/ + 1) + 2k{k + 1) ^^^^ 

for vector particles. Here /cq is the time component of the momentum and k^ is a momentum 
related to the ^-dimensional torus. 

The P2 term comes from the diagrams with two insertions. If the two time-momenta in 
these diagrams are denoted by p and p + g, wc first perform the integration over and then 
expand the result in q^. Using partial integration, these powers of can be transformed 
in time-derivatives acting on A and hence on B. There is also a diagram with an explicit 
dependence on P. It is the diagram with two insertions of the operator Wqi, one Qq and 
one Qi propagator. 
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Kn 


0566264741439181 


, — 


— 245345QQ8517Q565 




3 6686Q17Q8 14223 


K n 


n 50(170320300661(1 


K4 = 


-0.839359633413003 


K5 = 


-0.849965412245339 


K6 = 


-0.06550330854836428 


^^^ = 


-0.3617122159967145 




-2.295356861354712 



Table 1: The numerical values of the coefficients. 



Adding up the different contributions, we obtain the one-loop contribution to the kinetic 
term c^cf + dfdf in the lagrangian. Demanding the renormalised kinetic part to look just 
like the classical term gives us the finite part of the renormalisation (pBf): 
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g = I + 12;^ + 12^ ^"^(27^) + 

where kq can be found in table |1]. This renormalisation scheme can easily be related to 
other schemes like the MS scheme. 

The finite, renormalised effective potential becomes 

Veff=5-^cl(C'^) + ^eff' (40) 



with 



V^^icd) = V^^ic) + V^^id) + K,triX)triY) + KstTiXY), 
V^^{c) = /titr(X) + K2det{c) + Ksti^X) + K^tiiX'^) 
+K5 det(c) tr(X) + Kq tr^(X), 



(41) 



with the numerical values for Ki in table |l[ Note that the term in the effective potential 
along the tunnelling path uniquely determines the coefficient of the tr(X) det(c) term. The 
term can be obtained from combinations of the three independent invariants tr^(X), 
tr(X) tr(X^) and tr(X^). We choose to replace the term by tr^(X), which is the simplest 
of these from the viewpoint of the variational calculation. 



6 Variational Results and Conclusions 

We calculated the effect of the high-energy modes on the dynamics of the low-energy modes. 
It resulted in a renormalisation of the coupling constant and a correction to the potential 
in the effective hamiltonian. 

With the obtained hamiltonian, we can repeat the variational approximation of the 
spectrum. The results remain qualitatively the same: the lowest-lying scalar [j = 0) and 
tensor (j = 2) levels can be found in the same sectors as before, although the values of the 
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energy have changed. Results for the lowest glueball masses can be found in fig. 3. Around 
/ = 0.25 the mass ratio 1712+ /mQ+ is roughly 1.5, which compares nicely with the lattice 



results 



Just as for the lowest-order hamiltonian, we used Temple's inequality |TT[] to obtain 



lower bounds for the energy levels. This convinced us as before that we obtained accurate 
results. However, due to the more complicated structure of the potential, a larger number 
of basis vectors is required. 

The onset of the influence of the boundary can be seen at / = 0.2. One of the issues 
raised in P| was the level of localisation of the wave function around the sphaleron. This is 
related to the question whether the assumption is true that only the boundary conditions 
at and near the sphalerons are felt. We argued that this was determined by the rise of the 
potential in the transverse directions. The one-loop correction to the tr(y) term in the 
potential at the c-sphaleron, which can be expressed in ni, and ks, is such that it results 
in a lesser degree of localisation. Beyond / = 0.3 the approximation breaks down as can 
be seen for instance by the crossing of the scalar and tensor glueball. A fuller study into 
this localisation, as well as more details and results will be presented in the near future. 
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Figure captions 



Figure 1: Expansion of V^^\u) in the tunneling parameter u. We dropped 
the e and log(L) dependent parts. We have drawn the expansion up to order 
xi^, u^, and u^^. Longer dashes correspond to higher order expansions. 
The horizontal hne at u— 1 denotes the exact result at the sphaleron. 

Figure 2: Topology of the different Feynman diagrams. 

Figure 3: Glueball masses for = as a function of the coupling constant. 
The lower and upper drawn curves are the masses of resp. the first scalar (0^) 
and tensor (2"*") glueball. The dotted lines denote the perturbative results. 
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